Imaging and manipulating electrons in a ID quantum dot with Coulomb blockade microscopy 
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Motivated by recent experiments by the Westerveh group, which used a mobile tip to probe the electronic 
state of a segmented nanowire, we calculate shifts in Coulomb blockade peak positions, as a function of tip 
location, which we term "Coulomb blockade microscopy". We show that if the tip can be brought sufficiently 
close to the nanowire, one can distinguish a high density electronic liquid state from a Wigner crystal state by 
microscopy with a weak tip potential. In the opposite limit of a strongly negative tip potential, the potential 
depletes the electronic density under it and divides the quantum wire into two partitions. There the tip can push 
individual electrons from one partition to the other, and the Coulomb blockade micrograph can clearly track such 
transitions. We show that this phenomenon can be used to qualitatively estimate the relative importance of the 
electron interaction compared to one particle potential and kinetic energies. Finally, we propose that a weak tip 
Coulomb blockade micrograph focusing on the transition between electron number N = and A' = 1 states may 
be used to experimentally map the one-particle potential landscape produced by impurities and inhomogeneities. 
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I. INTRODUCTION 



Studies of nanoscale electronic structures hold important 
promise both as laboratories for few-body, interacting quan- 
tum mechanical systems and as technological testbeds for 
future classical or quantum computing technologies. Novel 
probe technologies ^ are very important for studying elec- 
tronic properties in nanoscale systems because they are often 
beyond the resolution of conventional imaging techniques like 
optical microscopy, and traditional transport measurements 
can only measure spatially averaged physical properties such 
as the conductance or the current. One scanning probe mi- 
croscopy (SPM)- technique utilizes a charged metallic tip to 
perturb the local electronic density in a nanoelectronic struc- 
ture while monitoring the resulting change in transport prop- 
erties. Using this technique one can obtain spatially resolved 
measurement of the electronic properties, including the local 
electron density and, in principle, the wavefunction itself in 
the case of a one-electron system2-(see discussions below). 
This imaging technique has been fruitfully applied to study 
the flow of ballistic electrons across a range of two dimen- 
sional heterostructures. 

Recently, a series of experiments^ applied SPM techniques 
to study quantum wires. In these experiments a segment of 
an InAs nanowire lying on top of a two-dimensional SiOv 
layer was isolated from the rest of the wire by two short 
InP segments, forming a one-dimensional quantum dot with 
lithographically defined boundaries. A negatively charged 
probe scanned controllably the two dimensional area around 
the wire and the conductance across the ID quantum dot was 
measured as a function of the probe location. Both the voltage 
of the probe and its height above the surface can also be inde- 
pendently varied. Motivated by these new experimental pos- 
sibilities, we turn to exact diagonalization techniques to study 
the conductance response of a few-electron quantum dot as a 
function of a spatially varied probe potential, in order to illus- 
trate the kind of information that can be extracted in the case 
of a system of several electrons. 




FIG. 1: Schematic geometry of Coulomb blockade microscopy of 
a quantum wire containing four electrons. When calculating the 
electron-electron interaction and the electron-tip interaction, we as- 
sume that the InP barriers have zero thickness, the InAs wire is in- 
finitely long, and the substrate layers extend to infinity in x and y 
directions. 



II. MODEL 

We consider a uniform InAs (dielectric constant e = 15.4, 
Bohr radius ag w 34nm) nanowire of radius R = lOnm, which 
lies in vacuum atop a SiOv (e = 3.9) layer lOOnm thick, sep- 
arating it from conducting doped bulk silicon (see Fig. [T]). 
Electrons are modelled as point charges traveling along the 
center axis of the wire, confined to interval — ^ < x < | 
by hard walls, representing the InP layers. We consider 
length L from llOnm to 500nm. The electron-electron in- 
teraction U{xi — X2) was calculated using the commercial 
finite-element program Comsol™to solve the classical Pois- 
son equation for a point charge on the axis of an infinite wire 
above in a substrate with the geometry described in Fig.[T] At 
short distance Ax, the potential was softened to account for 
the finite thickness of the electron wavefunction, by replac- 
ing Ax^' with [(Ax)^ +7?^]^'/^. Following the approximation 
used by Topinkai^, we model the negatively charged probe as a 
fixed point charge of strength ^ at a location ro relative to the 
center point of the wire. This gives rise to a one-body poten- 
tial V (x; ro , ^) for an electron on the wire axis at point x, which 
we again obtain by solving the Poisson equation (results are 
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FIG. 2: (Color Online) (Dotted)Interaction potential (/(Az). (Solid, 
long and short dashed) Tip V{x) potentials with a tip charge q = e 
and locations = (0,0, zq) where zq = 30nm,50nm, 100;im. 




FIG. 3: (Color Online) Electronic densities for a L = 1 lOnm wire and 
for a L = 500nm wire in absence of tip potential. Only the right half 
is shown, as the plot is symmetric about x = 0. We rescale the x and 
p{x) with wire length L to ease the comparison. 



shown in Fig.|2]l. 

In this paper, we diagonalize the exact ID many -body 
Hamiltonian with the Lanczos method^ for up to electron 
number N ~4: 

*2 N N ;-l 

2m f-'ipi 

(1) 

where ^ is the full many-body wavefunction, depending 
on the position Xi and spin a, of the electrons. To con- 
nect to the experimentally observable variables, we consider 
the Coulomb blockade peak positions of the transition from 
{N —l)to N electron ground states. The conductance through 
the quantum wire is maximum when the chemical potential 
difference between the lead and the wire, controlled by the 
voltage Vg on a back gate, is equal to the ground state en- 
ergy difference between the two states in question. We may 
write this condition as AE = En — En^\ = aVg + j3, where 
j8 is a constant and a is the proportionality constant between 
changes in the back gate voltage and the chemical potential in 
the quantum dot. We probe the electronic states in the quan- 
tum wire through the dependence of AE on the tip position 
ro and potential strength q. An interesting set of spatially re- 
solved information about the electrons in the wire can be ex- 
tracted from this function, and we call this method "Coulomb 
blockade microscopy". It is a special application of the "scan- 
ning probe microscopy" developed by the Westervelt group^. 
In calculations in this paper we focus on the transition from 
N = 3 to N = 4- electrons, but most of our conclusions are 
easily generalizable to other ground state transitions. 

Finally, we note that for four non-interacting electrons with 
spin in a wire of radius R = lOnm, when the dot length 
L > L„i = 18.2 nm, the lowest four single particle energy levels 
are all longitudinal modes. The shortest wire length we con- 
sider in this paper L = 1 10«m ^ L,„, we therefore expect the 
wires under consideration can be well approximated as strictly 
ID under the assumption of weak interaction effects. Indeed, 
the gap between the transverse ground state and first excited 



state for our wire is AE « l48t7ieV. The most confined ge- 
ometry we discuss is shown in Fig.|6] a L = 1 lOnm wire with 
four electrons under an extremely strong tips, squeezing them 
into the both ends of the wire. Even in that case, the total 
energy per electron is less than lAQmeV, not enough to cause 
an excited transverse mode, with full interaction effects taken 
into account. Thus we expect the quantum wire to be well ap- 
proximated by a ID model for all the parameters we explored 
in this paper. 

The existence of higher transverse modes will lead to a 
renormalization of the effective electron-electron interaction 
(generally a softening of the potential at short distances) even 
if the energies of the modes are well above the Fermi energy. 
This is because two electrons that come close together in the 
lowest transverse mode will mix, virtually, with states where 
the electrons are in excited modes. As our calculations use, 
already, a crude phenomenological cut-off at short distances, 
we do not include explicitly effects of this renormalization. 
Technically, the existence of higher transverse modes can also 
lead to three-body and four-body effective interactions, aris- 
ing from three-body and four-body collisions, but we do not 
expect such higher body terms to be important in the wires 
under consideration. 

We want to emphasize that although we choose a specific 
wire geometry and a specific form of interaction that model 
some characteristics of the setup of the Westervelt group's 
on-going experiments, most features of the Coulomb block- 
ade micrographs we discuss below are applicable to any ID 
quantum dot system under a mobile potential. Indeed, our 
discussions of the qualitative features of Coulomb blockade 
micrographs under both the weak-tip and the strong-tip lim- 
its rely only on the general properties of ID electronic sys- 
tems, independent of the specific geometry and interaction we 
adopt. Our numerical results mostly serve for illustrative pur- 
poses. 
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FIG. 4: (Color Online) Coulomb blockade micrographs for a ID 
dot with L = 500nm and tip charge q = 0.02e for three tip poten- 
tial shown in Fig|2] Again, the right half is shown. 



III. ELECTRONIC DENSITY 

In the absence of a probe potential V, both the N = 3 and 
N = 4 wires the electronic density profile p(x) undergoes a 
crossover as a function of L from a liquid state characterized 
by a Ikp Friedel oscillations to a quasi-Wigner crystal state 
characterized by a Akp density oscillation. Such a crossover 
from a liquid state to quasi-Wigner states with a decrease in 
density a very generic phenomenon for ID interacting fermion 
system^ *^. Indeed, for any interacting decaying no faster than 
x^^ at long distance a quasi-Wigner crystal state is known 
to emerge at low density^. This requirement for interaction 
will hold for a system with long screening length as compared 
to mean inter-particle distance, as is the case for our geom- 
etry when the screening doped silicon layer is relatively far 
(100/im) away. For our specific geometry and interaction, the 
crossover happens at around density p* « 35jLim^'. A Wigner 
crystallized density variation is shown in the dashed curve of 
Fig. [3]for a quantum dot of L = 500«m, whereas for L = 1 10 
the four electron density exhibits Friedel oscillations. 

IV. WEAK TIP LIMIT 

Now we introduce a weak tip potential, corresponding to a 
negatively charged tip of strength q = 0.02e, scanning above 
the center axis of the quantum wire along its direction (1,0,0), 
with the tip location vector ro = (xo,0,zo)- For a ID quan- 
tum dot of length L ~ 500«m, which as shown in Fig. |3]has 
Akp Wigner-crystal density variation, let us consider the three 
tip heights above the quantum wire, zo = 3Qnm,5Qnm, lOOnm, 
corresponding to the three tip potential shown in Fig. [2] The 
resulting Coulomb blockade peak position AE as a function of 
the tip coordinate xq along the wire, i.e. the Coulomb block- 
ade micrograph, is shown in Fig. 2] Clearly, in Fig. |2] the 
closer the tip approaches the wire, the more localized is the tip 
potential and a sharper tip potential make it easier to resolve 
the density variations, this is reflected in Coulomb blockade 
micrograph scans in Fig.|4] At zo = 30«m,50«TO from the tip 
to the center of the wire, the A-kp density oscillation of the 



FIG. 5: (Color Online) q = O.Ole weak tip limit for the yV = 3 to 
A' = 4 Coulomb blockade transition for densities shown in Fig. [3] 
The tip distance to the wire is zq = 30nm. Only the right half is 
shown. 
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FIG. 6: (Color Online) q = ie strong tip limit for the A' = 3 to W = 4 
Coulomb blockade transition for densities shown in Fig.[3]The tip to 
wire distance is 30nm. Only the right half is scanned. 



quasi-Wigner crystal state on the right can be detected in the 
Coulomb blockade micrograph, whereas when zo = lOOnm 
away, the tip potential becomes much too broad to resolve the 
fine features of the density oscillations. We note that although 
the resolution of the tip is largely determined by the distance 
Zo, the contrast of a Coulomb blockade micrograph, i.e. the 
magnitude of the 4kf variations in the micrographs, can be 
improved by modestly increasing the tip potential. 

By contrast, in Fig. [5] the L = llOnm micrograph at zo = 
3Qnm does not show features of Wigner crystal oscillations. 
However, this micrograph does not by itself give a clearcut 
indication of the absence of Wigner crystal order for L = 
llOnm. With the current interaction and tip parameters, one 
cannot observe the crossover from the Wigner crystal to the 
Friedel oscillations because it happens at a inter-particle spac- 
ing Ax « 30nm, below the resolution of the micrograph even 
at z = 30nm. To be more specific, we may define the onset 
of Wigner-crystal order for our four-electron system as the 
length L at which there first appears a local minimum of the 
mean density p{x) in the vicinity of x/L = 0.25. According 
to our calculations, this should occur at L=135 nm. However, 
with the tip at height 30nm, in the weak charge limit, the reso- 
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FIG. 7; Electronic densities in L = SOOnm, N = 4 quantum wire as a 
tip 9 = 8e scan througli < xq < 250«m. The density tliroughout the 
entire wire — 250«m < x < 250nm is shown. 



lution of the micrograph is of the order of 60nm, so we would 
not see a secondary minimum in the micrograph signal until 
L > 250«OT. 

To gain a more intuitive understanding of a weak tip 
Coulomb blockade micrograph, we observe that a weak tip 
only slightly disturbs the electron density as it scans across 
the wire, thus, a simple first order perturbation theory should 
be a good approximation to compute the ground state energy 
in the presence of the tip potential: 



E{r,q) -£(,{?, q)^ dxV{rQ;q,x)p{x), 



(2) 



where p{x) is the non-interacting ground state density and 
Eo{ro,q) is its energy. We have checked that for tip charges 
up to (7 = O.le the simple first order perturbation theory gives 
a decent fit to both the ground state energy and the Coulomb 
blockade micrograph. Since both the width and the center 
location of the tip potential V{ro;q,x) can be adjusted exper- 
imentally, the Coulomb blockade microscopy with a weak tip 
potential provides a flexible way to map the electronic densi- 
ties in a quantum dot. Indeed, as an example of such flexibil- 
ity, we find that one can improve the "contrast" of a Coulomb 
blockade micrograph; i.e., the prominence of the spatial vari- 
ations in micrographs like Fig. |4] as compared with the total 
energy shift £4 —£3, can be improved by slightly increasing 
the tip potential while still staying within the weak tip pertur- 
bative approximation. 



V. STRONG TIP LIMIT 

In the opposite limit of strong tip, the Coulomb blockade 
tip scans present a very different physical picture. In Fig. |6] 
we observe that irrespective of whether the electronic state is 
liquid or Wigner-crystal like as shown in Fig. [3] the Coulomb 
blockade micrographs show similar behavior: in the case of 

= 3 to = 4 transition, both the L = 1 lOnm and L = SOOnm 



FIG. 8: Schematic illustrations of partition of the four-electron state 
by the potential barrier under a strong tip potential. In (a) and (b), 
single-electron energy levels are shown for non-interacting electrons, 
for two positions of the tip. Because of spin degeneracy, we find ei- 
ther two electrons in each well, or all four in the same well, when 
the energy of the single-electron ground state on the right becomes 
higher than the first excited level on the left. In (c)-(e), we "incor- 
porate" interaction energy into "single electron levels" schematically 
by plotting the energy needed to add an additional electron. In draw- 
ing energy levels this way, it is clear that the spin degeneracy in the 
non-interacting case is lifted by electron interaction, and an addi- 
tional (3,1) partition will appear for the "energy level" arrangement 
in (d). Filled circles show electrons in occupied levels. 



wire show two relatively sharp peaks for a large tip charge 
^ = 8. This is in contrast with the case of a weak tip Fig. |5] 
where the Coulomb blockade micrographs show smooth spa- 
tial dependence as well as sensitivity to the electronic states in 
the absence of the tip potential. 

To understand the physics of this strong tip limit we note 
that the two sharp cusps in Fig. |6] represents discontinuous 
slope changes in the N = 4 electron ground state energy as a 
function of tip position xq. Similarly the deep valley in the fig- 
ure corresponds to a cusp in = 3 ground state energy. The 
origin of these three discontinuities in slopes can be seen in 
Fig- 111 In this limit, the negatively charged tip potential is so 
strong that it depletes the electronic density under it. Thus the 
tip creates an effective partition of the electrons in the wire 
into left and right sub-quantum dot. As shown in Fig. |71 as 
the tip move from the center to right of the wire with four 
electrons, the partitions of the electrons undergoes two abrupt 
transitions (2,2) — > (3, 1) — >• (4,0). These two transitions cor- 
respond to the two upward cusps shown in the N ~ 4 curve 
in Fig.|6l Similarly, the discontinuous slope change shown on 
the N = 3 curve of the Fig.|6]corresponds to the transition be- 
tween the (2, 1) — > (3,0) partition of the ground state. Thus 
the three discontinuities seen in the Coulomb blockade mi- 
crographs in Fig. |6]coiTespond to, alternately, the transitions 
between the integer partitioning of total electron numbers in 
the N ^3 and N ^ 4 system. The upward slope of the curve 
near . to = reflects an additional downward cusp at the origin, 
due to the transition (1,2)— >(2,1) in the N ^3 wire. 

To better understand the transitions between different parti- 
tions, let us consider the transitions in a model of = 4 elec- 
trons, with spin, which have no Coulomb repulsion between 
them but interact with a repulsive tip potential. The scenario 
is illustrated in Fig. [8] panels (a) — (b). When tip is near the 
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FIG. 9: (Color Online), q = 0.6e and q= 1.8e intermediate tip po- 
tentials for W = 3 to A' = 4 Coulomb blockade transition in a wire of 
length L = 180wm. The tip distance to the wire is zq = 30nin. Only 
the right half is shown. 



FIG. 10: (Color Online) Electronic density for a L = IWnin wire 
under three tip potentials, with effective charges q = Oe, q = 0.6e 
and q = 1.8e. Only the right half is shown, as the plot is symmetric 
about X = 0. We rescale the x and p(x) with the wire length L. 



center of the wire, the electrons are partitioned (2,2) and both 
electrons in each side reside in the single-particle ground state. 
As the tip moves rightward, the energy levels rise in the right 
partition and fall in the left. When the first excited level on the 
left partition crosses the ground state on the right, both elec- 
trons will move to the left partition. Therefore, contrary to the 
interacting case, there is no energetically favorable state of 
(3,1) partitioning in the non-interacting system. In the N = 3 
case, the (2, 1) partition is not affected by this, and for non- 
interacting system the transition (2, 1) — > (3,0) will coincide 
with the transition in = 4, so the Coulomb blockade micro- 
graph will show only a single peak. This analysis can also 
be generalized to a wire containing multiple non-interacting 
electrons, such that all the {odd, odd) partitionings of elec- 
tron number will be missing. 

As schematically illustrated in Fig. [8] panels (c) — (e), 
when we take electron interaction into consideration, the one 
and two electron state would no longer be degenerate in either 
partition, so contrary to the non-interacting scenarios in panel 
(a) — (b), here a (3,1) partition can survive as an interme- 
diate stage between the (2,2) and (4,0) partitions. With the 
non-interacting case in mind, we postulate that the distance 
between two peaks in = 3 to = 4 micrographs, coitc- 
sponding to the tip positions where (3,1) partitioning in the 
N ^ 4- wire is stable, can serve as an indicator of the rela- 
tive importance of the interaction energy versus the sum of 
kinetic and single particle potential energies. The less impor- 
tant interaction is compared to single particle energies, the less 
splitting would the one and two particle energies be, and the 
smaller is the region of stable (3,1) partition. This can be seen 
in Fig. |6] The potential energy should have a larger share in 
the total energy in the longer wire with lower electronic den- 
sity, and indeed we observe that the longer wire has a wider 
distance between the two peaks marking (2,2) — (3,1) and 
(3, 1) (4,0) transitions. 



VI. INTERMEDIATE TIP CHARGE 

We have carried out calculations with various tip charges 
intermediate between the two limits discussed above. In gen- 
eral, the larger the charge on the tip, the more readily one 
sees the secondary minima in the micrograph signal, which 
are seen in Fig. |6]for a charge q ~ 8e, even at L = llQnm. 
As one illustration, for a wire of length 180 nm, with a tip 
height of 30 nm, we find that the micrograph signal shows a 
secondary minimum at a:/L « 0. 15 for q = 1 .8e, but shows no 
secondary minimum when q = 0.6e, as seen in Fig. |9] How- 
ever, if we calculate the electron density in the wire when 
the tip is over the center of the wire (xq = 0), we find that 
a tip charge of ^ = 0.6e is enough to substantially modify the 
density relative to the density in the absence of the tip. As 
seen in Fig. [TOl the electron density below the tip, at x = 
is reduced by a factor of three relative to the density with no 
tip charge. Nevertheless, the oscillations seen in p{x) remain 
qualitatively similar to structure seen in the absence of the tip. 
For example, the charged tip only pushes out the position of 
the secondary minimum in the density from x/L = 0.25 to 
x/L = 0.3. 

In general, when we increase the tip potential to intermedi- 
ate values, the resolution of our micrographs improves com- 
pared with the weak tip limit. This is evident in the fact that at 
q= 1 .8e we can already see signs of Wigner crystallization at 
L = 180«m, in contrast with the weak tip case with q = 0.02e 
where we can only detect quasi-Wigner crystal at L = 250«m. 
On the other hand, atq= 1 .8e there is no signature of Wigner 
crystallization when the electrons are in a liquid state in a 
L — llOnm wire. Thus we do not have "false positive" sig- 
nature of quasi-Wigner crystal, in contrast with the strong tip 
limit described in Sectior(y] where the micrographs show A^ 
peaks for a A^-electron wire regardless whether, in the wire in 
the absence of the tip, the electrons are in a liquid or a quasi- 
Wigner crystal state. 

Because of these two characteristics, an intermediate tip po- 
tential may help an experimentalist to reliably detect the pres- 
ence of a quasi-Wigner crystal state in a shorter, higher density 
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wire closer to the crossover from a liquid state. 

VII. SINGLE ELECTRON 

Beyond the system of interacting electrons discussed above, 
a possible further application of Coulomb blockade mi- 
croscopy is to experimentally "map" the rugged potential 
landscape produced by wire inhomogeneities and charged im- 
purities in the substrate. One would focus on the transition 
from N ^Oto N = 1 state, in which case the Coulomb block- 
ade micrograph would reveal information about the single par- 
ticle density. By inverting the transformation in Eq. |2] one 
may be able to approximately obtain the single particle ground 
state density p{x). In the absence of an external magnetic 
field, the ground state wavefunction \lf{x) has no nodes and 
can be chosen to be yf{x) = (p(x))'''^ It is then straightfor- 
ward to invert the Schrodinger's equation to extract the poten- 
tial landscape from the single particle wavefunction. 



the other, and the accompanying Coulomb blockade micro- 
graph can indicate the transitions between different partition- 
ings. Furthermore, a feature of the resulting micrograph, the 
distance between peaks marking the {odd, odd) partitioning, 
can serve as an indicator of the relative strength of the inter- 
action. 

In this paper we have chosen extreme values of the tip 
charge q to illustrate the physics in the two limits. However, 
our calculations show that the discussions above hold true for 
a wider range of moderately small and large values of q. 

To obtain a quantitative description of the energy shifts ex- 
pected in Coulomb blockade microscopy, particularly in the 
intermediate coupling regime, we see that it is necessary to 
perform a realistic calculation, which takes into account both 
the electron-electron interaction and the non-linear effects of 
the charged tip on the electronic state of the wire. If one is 
prepared to carry out such a calculation, however. Coulomb 
blockade microscopy can be a powerful probe of interaction 
effects in the wire. 



Vni. SUMMARY 

In summary, in this paper we show that tracking the peak 
position shift as a charged mobile tip moves above and across 
a nanowire, a technique we term Coulomb blockade mi- 
croscopy, can reveal spatially-resolved information about the 
electronic density and states of a quantum ID dot. A weak tip 
potential can serve as a probe with a tunable width, to reveal 
the spatial distribution of the electronic density in the wire. 
A strong tip potential that depletes part of the wire can be 
used to manipulate individual electrons from one partition to 



IX. ACKNOWLEDGEMENTS 

We would like to thank Erin Boyd, Halvar Trodahl and 
Jesse Berezovsky and especially Bob Westervelt for helpful 
discussions. This work is supported in part by NSF grants 
PHY-0646094 and DMR-0906475. Numerical work was per- 
formed in part at the Center for Nanoscale Systems, a mem- 
ber of the National Nanotechnology Infrastructure Network 
(NNIN) supported by NSF award ECS-0335765. JQ is also 
supported in part by NIM and DFG through SFB 63 1 . 



^ O. M. Auslaender, A. Yacoby, R. de Picciotto, K. W. Baldwin, 
L. N. Pfeiffer and K. W. West Science, 295(5556):825-828, 2002. 

2 M. A. Topinka, R. M. Westervelt and E. J. Heller Physics Today, 
56(47), 2003 and refs. therein. 

^ P. Fallahi, A. C. Bleszynski-Jayich, R. M. Westervelt, J. Huang, 
J. D. Walls, E. J. Heller, M. Hanson and A. C. Gossard, Proc. 27th 
Int. Conf. on Phys. of Semiconductor, AIP, 2005, p. 779. 

^ A. C. Bleszynski-Jayich, L. E. Froberg, M. T. Bjork, H. J. Trodahl, 
Lars Samuelson and R. M. Westervelt, Phys. Rev. B, 77, 245327 



2008. 

^ M. A. Topinka. PhD thesis, Harvard University, 2002. 

^ J. K. CuUum and R. A. Willoughby. Lanczos Method for Large 

Symmetric Eigenvalue Computation, Vol. 1. Birkhauser Boston 

Inc., Boston, MA, 1985. 
^ G. A. Fiete, J. Qian, Y. Tserkovnyak, and B. I. Halperin Phys. Rev. 

S, 72, 045315, 2005. 
^ K. A. Matveev, Phys. Rev B, 70, 245319, 2004. 



